Home > Publications database > Numerical integration of the reactor kinetics equations with exponentially fitted implicit methods |
Book/Report | FZJ-2018-00881 |
1976
Kernforschungsanlage Jülich, Verlag
Jülich
Please use a persistent id in citations: http://hdl.handle.net/2128/16869
Report No.: Juel-1307
Abstract: This report treats the numerical integration of the reactor kineties equations for one and two prompt neutron groups by means of advanced numerical techniques being able to solve extrernely stiff initial value problems accurately at lange discretization intervals. The numerical algorithms being considered belong to a dass of exponentially fitted irnplicit methods with unlimited stability fange. Fitting methods being improved compared to existirtg approaches male an accuracy with lange discretization intervals possible, which, at comparable computational effort, may be several orders of magnitude baffer than the accuracy of established methods which are designed for the treatment of stiff problems. This is demonstrated frone results of numerical experiments which compare the proposed rnethods with known A-stable and strongly A-stable irnplicit methods without exponential fitting. The improved methods require an adaptation of the algorithm to the properties of the system of differential equation. The general considerations which have to be kept in view are explained so far as necessary for transmitting the technique to other stiff systems.
The record appears in these collections: |